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ABSTRACT 



We have compiled a sample of 16 orbits for very low-mass stellar (< 0.1 Mq) and 
brown dwarf binaries, including updated orbits for HD 130948BC and LP 415-20AB. 
This sample enables the first comprehensive study of the eccentricity distribution for 
such objects. We find that very low-mass binaries span a broad range of eccentricities 
from near-circular to highly eccentric (e ~ 0.8), with a median eccentricity of 0.34. We 
have examined potential observational biases in this sample, and for visual binaries we 
show through Monte Carlo simulations that if we choose appropriate selection criteria 
then all eccentricities are equally represented (< 5% difference between input and output 
eccentricity distributions). The orbits of this sample of very low-mass binaries show 
some significant differences from their solar-type counterparts. They lack a correlation 
between orbital period and eccentricity and display a much higher fraction of near- 
circular orbits (e < 0.1) than solar-type stars, which together may suggest a different 
formation mechanism or dynamical history for these two popu l atioii s. Very low-mass 
binaries also do not follow the distribution of lAmbartsumianI ( 19371 ). which would be 
expected if their orbits were distributed in phase space according to a function of energy 
alone (e.g., the Boltzmann distribution). We find that current numerical simulations of 
very low-mass star formation do not completely re produce the observed properties of 
our binary sample. The cluster formation model of iBatd (j2009al ) agrees very well with 
the overall e distribution, but the lack of any high-e (> 0.6) binaries at orbital periods 
comparable to our sample suggests that tidal damping due to gas disks may play too 
large of a role in the simulations. I n con trast, the circumstellar disk fragmentation 
model of IStamatellos &: WhitworthI (120091 ) predicts only high-e binaries and thus is 
highly inconsistent with our sample. These discrepancies could be explained if multiple 
formation processes are responsible for producing the field population. 
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Introduction 



Binary systems have long been used as unique probes into the origins of stars, as the dy- 
namical imprint of their formation and subsequent evolution are recorded in their orbits. In fact, 
decades before it was known that nucleosynthesis powers stars, well-dete rmined binar y orbits were 
commonplace (e.g., > 200 orbits were published in the compilation of lAitkenl Il918l ). Thus, the 
study of binary orbits has a rich heritage, predating much of modern astrophysics. However, only 
very recently has it been possible to measure orbits for the lowest mass products of star formation: 
stars at the bottom of the main-sequence (< 0.1 Mq) and brown dwarfs. The sample of orbits 
for such objects has been growing rapidly in recent years as astrometric monitoring has begun to 
yield orbits for the shortest period visual binaries know n (P ^ 10 years), which were discovered 
by high-resolution in i aging surveys 5-10 year s ago (e.g., iReid et al.ll200ll . l2006l : IBouv et al.l 120031 : 
Burgasser et"al]l2003l . bood : Iciose et ahlbood ). 



Some properties of very low-mass binaries, such as orbital semimajor axes, mass ratios, and 
multiplicity fractions, can be studied in a s tatistical fashion using a large sample without any knowl- 
edge of orbital motion (e.g., lAllenI 120071 : iBurgasser et al.M2007l ). However, only well-determined 
orbits can provide a measure of eccentricity (e), which is directly related to the dynamics of the 
formation process. Because of the long history of stellar orbit studies, many authors have weighed 
in on the subject of the eccentricity distribution^ The most recent comprehens ive work on the 
eccentricities of solar- type field binaries was done b y IPuauenriov fc MavoJ (1991), and their pre- 
main-sequence counterparts have been reviewed by Mathieu ( 19941 ). Now for the first time it is 



possible to bring similar observations to bear on the extreme low-mass end of the binary popu 
latior i. Some of the proposed formation mechanisms for such low- mass qbiect s , such as ejection 



rope 
3E 



(e.g.,lReipurth &: Clarkdl200ll ) and gravitational instability (e.g.. lRice et al. 



2003 



Stamatellos et al 



20071 ). may leave a distinct signature in the orbits of binaries. Here we present the first compre- 
hensive look at what orbital eccentricity measurements reveal about the nature of very low-mass 
stars and brown dwarfs. 



2. Updated Orbit Determinations 



We first updated published orbit determinations for very low-mass binaries that have new data 
before compiling our sample for this study. For two visual binaries (HD 130948BC and LP 415- 
20AB), we present new astrometry from our ongoing Keck adaptive optics (AO) monitoring program 
in order to refine their orbital elements. Our procedure for obtaining, reducing, and analyzing our 
natural guide star a nd laser guide star AO images is described in detail in our previous work (e.g., 
Dupuy et al. 2009bl l3) . and the binary parameters derived from our new observations are shown in 



^According to lAmbartsumianl (1937|), "the distribution of eccentricities of double star orbits has been discussed 
often enough." 
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Table m We used the Markov Chain Monte Carlo (MCMC) approach described bv lLiu et al.1 (|2008l ) 
to determine the posterior probability distributions of the orbital elements for these binaries, and 
in Table [2] we give their median and 68.3% co nfidence limits . For H P 130948BC, we combined our 



new astrometry with the epochs published bv iDupuv et al.l (l2009bl). We update d this astrometry 
to account for the newer astrometric calibration of NIRC2 bv iGhez et al.l (j2008l ). which amounts 
to a +0?26 offset to our previously published position angles. The extended orbital phase coverage 
results in a factor of w4 improvement in precision compared to the published orbital parameters. 

Some binaries have multiple orbit determinations in the literature, and in these cases we 
used our published results in the following analysis. In most cases published orbits are in good 
agreement with each other, but a few binaries warrant special i nention. In the case of LP 415-20AB, 
our astrometry results in a very different orbit than found bv lKonopackv et al.l (120101) . As a test, 
we fitted their three Keck epochs with the published Gemini epoch from Isiegler et al.l (j2003l ) and 
found an orbit consistent with the one they reported (e = 0.9 it 0.1). However, after incl uding our 
four Keck epochs we found that the 2007 December astrometry of lKonopacky et al.l (|2010l ) is highly 
discrepant with the ensemble of measurements, resulting in a x'^ of 21.5 for 9 degrees of freedom 
(DOF). If we exclude this one measurement, we find good agreement between the data sets, and 
the best-fit orbit has a reasonable \^ (6.83 for 7 DOF) and significantly different orbital parameters 
(e.g., our newly derived eccentricity of 0.708j]QQjg is 1.9( 7 discrepant with their 0. 9ifc0.1). It appears 



that with only four epochs (i.e., 1 DOF) the orbit fit of iKonopackv et al.l (|2010j) was vulnerable to 
a single discrepant astrometric measurement leadin g to an incorrect orb it solution. We also note 
that our pubHshed orbi t for 2MASS J220 6 -204 7AB (jPupuy et al.ll2009al ) is nominally inconsistent 
with that published by lKonopacky et al.l (|2010l ). with all orbital parameters discrepant by ~2-4(T. 
Using only their data, we have found a lower orbit solution than they report, and although it is 
more consistent with our published orbit, it still has an unrealistically large (30.3 for 9 DOF). 
Thus, we choose to use our pubhshed orbit for 2MASS J2206-2047AB (x^ = 7. 49 for 7 DOF) in 



the fo llowing analysis. Finally, the orbit for 2MASS J 1534— 296 2 AB publi s hed b v lKonopackv et al 



(j2010l ) is marginally inconsistent with the original orbit from iLiu et al.l (j2008l ) . With the latest 
astrometry the period and semimajor axis are smaller by l.Scr, while the ec c entric ity and mass are 
consistent within \a. We use the more recent orbit from IKonopackv et al.l (|2010l ) in the following 
analysis. 



3. Very Low-Mass Binary Sample 

We have constructed a sample of very low-mass (Mtot < 0.2 A/q) binary orbits, including both 
visual and spectroscopic binaries. Because it is imperative for our purposes that this sample be 
as unbiased as possible with respect to eccentricity, in the following analysis we consider potential 
sources of selection bias, the nature of which varies between different types of binaries. Below we 
describe how we have selected our orbit sample as well as the tests we have performed to quantify 
the level of eccentricity bias present in this sample. 
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3.1. Visual Binaries 



In Table [3l we present a list of all 20 visual binaries that have published orbit determinations 
and a total system mass < 0.2 Mq. For each binary we list its distance, semimajor axis (a), 
and eccentricity (e), as these are the physical parameters that, along with three viewing angle 
parameters, determine the size and shape of a binary's orbit on the sky. Thus, these are the three 
key parameters that must be considered when characterizing selection biases. In most cases we 
have used orbit parameters as published and distances determined from trigonometric parallax 
measurements. When parallax measurements were not available, we estimated distances based 
on the primary components' resolved photometry and spectra l types us i ng th e updated spectral 
type-absolute magnitude relations described in Appendix B of iLiu et al.l (|20ld ). 



For all orbits published by other authors, we have created our own Markov chains using the 
published astrometry and the same analysis method described in Section [2l However, only in 
the cases where the astrometric orbit fits yield a reasonable can we trust the MCMC analysis 
to provide accurate posterior distributions of the orbital parameters. In cases where is un- 
reasonably large, the parameter space is not fully explored and uncertainties can be significantly 
underestimated by the MCMC method (or in fact by any fitting m ethod). We find such la rge 
values (reduced = 2.8-120) for five of the orbits published bv iKonopackv et al.l (j20ld ). For 
three of these orbits they combined astrometry with radial velocities (2MASS J1847+5522AB, 
2MASS J1750+4424AB, 2MASS J1426+1557AB), which may have enabled them to achieve or- 
bit fits with smaller than we found using the astrometry alone. For the other two orbits 
(2MASS J0850+1057AB and 2MASS J1728+3948AB), it is unclear why we are unable to repro- 
duce their lower x^ values from th e same input data (astro metry alone) . For all five of these orbits 
we use the parameters derived by lKonopacky et al.l (|2010l ). (Note that for reasons unrelated these 
discrepancies, none of these five binaries ultimately passes the selection criteria that we define below 
for our unbiased eccentricity sample.) 



3.1.1. Selecting a Minimally Biased Sample 



With a complete set of the known very low- mass visual binary orbits in hand, we now consider 
how selection biases may affect the eccentricity distribution of this sample. Visual binaries must 
pass three independent selection criteria before having their orbits determined: (1) they must be 
initially resolved by a high-resolution imaging survey; (2) they must be included in subsequent 
astrometric monitoring programs; and (3) their observed relative motion must yield a unique orbit 
fit. The first selection cut can result in biases on orbital parameters, as has been shown for direct 



imaging searches for extrasolar planets (e.g., see [Brown 



20041 ) . The second cut results from the 



fact thi .t . impractical to monitor the orbits of all .100 Known very low-m..= binarre.i Rather. 



^See www.vlmbinaries.org, last updated 2009 July 28. 
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astrometric monitoring programs are inevitably optimized for measuring as many orbits as possible 
with limited observing resources over a limited amount of time by focusing on the binaries likely to 
have the shortest orbital periods. Thus, the selection criteria used to define the monitoring sample 
can result in strong biases on the orbit parameters. Because the orbital period must be estimated 
from the projected separation of the binary, binaries with wide separations are much less desirable 
to monitor. This introduces an eccentricity bias because more eccentric binaries at a given period 
and mass will, on average, have wider projected separations. We will therefore attempt to define a 
subset of visual binaries that has minimal eccentricity bias due to these selection criteria. 

• Discovery bias. We first consider the effect of "discovery bias" which causes the population 
of binaries resolved by imaging surveys to deviate from the true underlying eccentricity distribution 
of all binaries. For a given imaging surveys, the extent of this bias depends on how the semimajor 
axis (a) distribution compares to that survey's resolution limit, or inner working angle (IWA). For 
example, if the a distribution peaks inside the IWA, the discovery of higher eccentricity binaries 
is favored, as such binaries spend more time at wider separations compared to more circular or- 
bits. To quantify this bias, we have performed a simple Monte Carlo simulation in which a was 
held fixed, while eccentricity was distributed uniformly from to 1, viewing angles were drawn 
from appropriate random distributions {p{i) oc sini; uniform in cj and fJ), and the single-epoch 
observations were distributed uniformly in orbital phase. Binaries were considered detected if their 
projected separation at the time of observation was above the resolution limit. The impact on the 
uniform input eccentricity distribution is illustrated in Figure [TJ When the resolution limit falls at 
or outside the fixed semimajor axis (IWA > a), high-e orbits are strongly preferred as they spend 
the most time wider than the resolution limit. When the resolution limit falls inside the fixed 
value of a (IWA < a) both highly eccentric and circular orbits are preferred, with intermediate 
eccentricities being selected against. This is because in these cases any sufficiently eccentric orbit 
can be made unresolvable by an unfavorable viewing angle, but this can be compensated by the 
fact that the most eccentric orbits spend extended periods of time at very wide separations. Thus, 
high-e orbits are more likely to be detected than intermediate-e orbits, as they spend the most time 
widely separated. 

To predict the exact form of the eccentricity discovery bias for our sample would require an 
assumption for the underlying semimajor axis distribution as well as careful accounting of the 
resolution limits of the numerous imaging surveys that have contributed to the current sample of 
very low-mass binaries. However, for our purposes we simply want to ensure that our sample of 
orbits has an acceptably small amount of bias. Thus, we need only chose a cutoff in the ratio 
of IWA /a that is small enough to minimize the bias while including as many binaries with orbit 
determinations as possible. In Table [3l we give our estimate for the IWA of the imaging data used 
to discover each visual binary as well as the ratio IWA/a. We have chosen a cutoff of IWA/a < 0.75, 
which Figured] shows produces a relatively unbiased selection over a full range of eccentricities from 
to 1. This cutoff excludes three binaries: LP 415-20AB (IWA/a 1.28), 2MASS J0920+3517AB 
(IWA/a w 0.91), and 2MASS J2140+1625AB (IWA/a w 0.83). The highest remaining IWA/a 
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value in our sample is 0.68 for LP 349-25AB. While it is unfortunate to have to exclude any orbits 
from our sample, we emphasize that this is a very important step in keeping potentially large 
biases from impacting the following analysis. The case of LP 415-20AB is a prime example of the 
problem with discovery bias on eccentricity. This binary was disco vered at a proiected separation 
of 119 ±8 mas, right at the resolution limit of the imaging survey of lSiegler et al.l (j2003l ). Its actual 
semimajor axis turns out to be 93.5^2 3 mas, so it could have only been discovered by that survey 
if it had an eccentric orbit, as it indeed does (e = 0.708^qqJ4). Thus, while LP 415-20AB does 
tell us that very low-mass binaries can be quite eccentric, it cannot convey any other information 
about the eccentricity distribution since the only reason it was discovered is because of its large 
eccentricity. 

• Monitoring sample bias. Next, we consider how to prevent the sample selection of 
orbital monitoring programs from biasing the resulting eccentricity distribution. This is essentially 
a matter of choosing a range of semimajor axes for which the monitoring sample is complete at all 
eccentricities. The lower limit in a is simply defined by the visua l binary orbit with the smallest 
semimajor axis; this is Gl 569Bab with a = 0.923 it 0.018 AU (jPupuy et al.ll20ld ). There are 
also no known visual binaries with projected separations (p) at discovery smaller than 0.9 AU, so 
we define this as the lower limit of our sample. The upper limit should be chosen to accurately 
reflect the completeness of orbit monitoring surveys that tend to focus on the smallest separation 
binaries, while also including as many of the resulting orbit measurements as possible. For the many 
binaries without orbit determinations, we must use their minimum possible semimajor axis, which 
would be the case of a highly eccentric, face-on orbit caught at apoastron (physical separation of 
a(l -|- e)). Since bound orbits have < e < 1, the minimum possible semimajor axis for any binary 
is therefore half its projected separation (p/2). Thus, some binaries that seem to be very wide (i.e., 
poor monitoring targets) could in fact have much smaller semimajor axes but just be very eccentric. 
When choosing an upper limit in a, we must therefore ensure that orbit monitoring programs are 
complete up to projected separations of 2a, otherwise we would introduce a strong selection bias 
against high eccentricities. W e choose an upper lim it in a of 3.7 AU in order to includ e the binaries 



LHS 1070BC (3.57±0.07 AU: ISeifahrt et al.ll2008l ) and LHS 1901AB (3.70±0.16 AU: IDuduv et al 



2OIOI). The next widest binary that would also pass our IWA/a cut is 2MASS 0850-hl057AB 



(4.8;?:^ AU, e = 0.64 ± 0.26). Our chosen 3.7 AU upper limit means that binaries up to projected 
separations of 7.4 AU should not have been intentionally excluded from orbit monitoring programs. 
To the best of our knowledge this is the case, as at least a second epoch has been obtained for 
all binaries with < 7 .4 AU, either by our own Keck program or others (e.g 



Bouy et al. 



2008 



Konopacky et al.ll2010l ). A slightly larger upper limit of 4.0 AU would require several more binaries 
to have undergone orbit monitoring that have not, while adding no additional binaries with orbit 
determinations to our sample. 

• Monte Carlo simulations. With these selection criteria of 0.9 < a < 3.7 AU, p < 7.4 AU, 
and IWA/a < 0.75 in hand, we have performed Monte Carlo simulations to assess the level of 
eccentricity bias that is expected for our sample. For this range of semimajor axes, the very low- 
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mass binary parameter distributions derived by lAllenI (|2007l ) are appropriate. The semimajor axis 
distribution is lognormal, with a peak at log(a/AU) = 0.86 and uioga = 0.28. To convert a from AU 
to arcseconds, we randomly drew distances in a Monte Carlo fashion incorporating the selection 
effects of the seeing-limited surveys that originally discovered the binaries as unresolved sources. 
We began with objects distributed uniformly in volume. W e assigned magnit udes to these objects 
according to the empirical ^C^-band luminosity function of I Cruz et al.l (120071 ) which is valid from 
9.5 < Mj^^ < 13.0 mag (i.e., ^M8 to L8, encompassing the spectral type range o f the bulk o f 
our sample) . We further added flux due to unresolved companions according to the lAUenl (j2007l ) 



binary frequency of 20% and mass ratio distribution of q^'^; we converted m ass ratio to flux ratio 



using the L oc M^'^ theoretical power-law relation of iBurrows et al.l ( 200ll ). The final distance 
distribution was then determined by magnitude cuts on this simulated population. We tested both 
strong (Ks < 12. ma g) and weak (Ks < 15.0 mag) magnitude cuts equivalent to those used by 



Gizis et al 



((20001) and lKirkpatrick et all (|2000l ) in mining 2MASS for late-M dwarfs and L dwarfs, 
respectively. In the latter cas e, we applied an addition al Mj^^ > 11.0 mag cut to approximate the 
J — K > 1.7 mag color cut of iKirkpatrick et al.l (j2000l ) that excluded brighter, bluer M dwarfs; we 



based this cut on the color-absolute magnitude diagrams of iDahn et al.l (120021 ). 



The strong magnitude cut for the brighter late-M sample resulted in a distance distribution 
peaking at ~20 pc, while the weaker cut for the fainter L dwarf sample resulted in similarly shaped 
distribution peaked at f«40 pc. We then randomly generated an observed projected separation 
(p) for each binary in a similar manner as described earlier (i.e., random viewing angles and 
observing times, and a uniform eccentricity distribution). Figure [T] shows the resulting eccentricity 
distribution after applying the selection criteria of 0.9 < a < 3.7 AU, IWA/a < 0.75, p < 7.4 AU, 
and p > IWA (i.e., the binary was resolved). For both samples we conservatively assumed the 



20031). and for the L dwarfs this is 0^.^07 from both HST and Keck surveys 



20od; 


Sieeler et al. 


, Bouv et al. 


2003; 



Reid et al.ll2006l : ISiegler et al.ll2007l ). By assuming the worst-case IWA we find the largest possible 
amplitude of the eccentricity bias, since smaller values of the IWA will always yield less bias (i.e., 
more binaries will always be uncovered for a smaller IWA). The resulting maximal bias we find 
from our Monte Carlo simulations has an amplitude of less than ±5% for both the late-M and L 
samples. As expected, the level of eccentricity bias introduced by our selection criteria is negligible, 
since they were chosen to minimize the eccentricity bias, and as shown in Figured] there is no strong 
systematic trend in the bias. Nearly circular (e < 0.2) and highly eccentric orbits (e > 0.8) are 
slightly preferred (by ~4%) as compared to intermediate eccentricities. We neglect this bias in the 
following analysis. 

• Orbit fitting bias. Finally, we consider one additional selection effect present in our 
sample: each binary must have undergone astrometric motion suitable for yielding a robust orbit 
determination. This can potentially introduce an eccentricity bias if orbit fits are more easily 
determined for certain ranges of eccentricity. For example, it is easy to imagine that an extremely 
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eccentric binary would present a challenge to astrometric monitoring as it spends most of its time 
at a fixed projected separation before moving rapidly inward in nearly a straight line. Such an orbit 
would be difficult to observe at just the rig ht times needed for a good fit and would also appear to 



be less appealing for continued monitoring. [Harrington &: MiranianI (|l977l ) investigated this bias by 



generating 10^ random orbits and classifying the expected quality of the orbit fit by eye, relying on 
their extensive experience working with such astrometric data. Such an ad hoc method is necessary 
because in the marginal cases of visual binary orbit fi tting the process of ininimiz ation is highly 



nonlinear and strongly dependent on the initial guess. [Harrington Sz MiranianI (| 19771 ) found that on 
average orbits with e < 0.6 all had the same likelihood of yielding a good fit, so they assigned this 
eccentricity range a normalized completeness factor of 100%. They found that this completeness 
drops to 77% for eccentricities of 0.6-0.8 and 42% for e > 0.8. In our sample, we have one visual 
binary with e > 0.8 and none in the 0.6-0.8 range. In the following analysis we consider the impact 
of applying correction factors of 1.3x and 2.4x due to this orbit fitting bias. 

• Final visual binary sample. In Table HI we list the orbital periods, eccentricities, and 
masses of the eleven visual binary orbits from Table [3] that meet our eccentricity bias- minimizing 
criteria of having 0.9 < a < 3.7 AU and a discovery TWA/ a < 0.75. (Note that we list orbital 
period here rather than semimajor axis in order to compare these orbits directly to the unresolved 
spectroscopic binary orbits.) Although this nominally includes only a little over half of the 20 
published visual binary orbit determinations, in fact the majority of the remaining orbits would not 
contribute much information as they are not as well constrained (see the orbit quality metrics given 
in Table [3]). For example, all orbits in our bias- minimized sample happen to have la eccentricity 
confidence limits of ^0.1 (median 0.005). However, of the nine orbits not included in the sample, 
six of them have much more poorly constrained orbits {la confidence limits of 0.1-0.5), with a 
median error of 0.26 that is > 50x larger than the median for the other orbits. In fact, the reason 
why these six orbits are poorly determined is the same as why they were excluded by our bias 
minimizing criteria: they are all the widest (i.e., longest period) binaries, and this results in them 
having orbit monitoring observations that cover much less orbital phase (as little as 0.2% and 1.3% 
of the orbit is covered in the two most extreme cases). 



3.2. Spectroscopic Binaries 



To date, five unresolved very low-mass binaries have published orbit determinations. Two 
of these were discovered via radial velo city (RV) monitoring as single-lined spectroscopic binaries 



(SBls): Cha Ha 8 (e = 0.59 ± 0.22: I.Toergens fc Miilleil 120071'). and 2MASS J0320-0446 (e 



0.067±0.015; Blake et al. 2008 ). Shen Turner ( 20081 ) have examined the selection bias introduced 



by RV surveys on the eccentricity distribution through Monte Carlo simulations. They found that 
as long as the RV semiamplitude (-f^i) of an SBl is ^3x larger than the rms of the residuals (a), 
its eccentricity does not greatly impact its likelihood of discovery. As shown in their Figure 4, 
SBls have a constant detection efficiency up to e = 0.6, beyond which there is a slight downturn. 
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However, even at the highest eccentricities (0.85 < e < 1.0) the detection efficiency is only f«20% 
less than for e < 0.6 orbits. The reason for this lack of a strong eccentricity bias in RV data is 
because eccentricity acts to boost the RV semiamplitude {Ki oc l/Vl — e^), counteracting the fact 
that eccentric orbits are more difficult to adequately sample with times series data. The two SBls 
in our sample, Cha Ha 8 and 2MASS J0320— 0446, have large semiamplitudes relative to their rms 
{Ki/a = 11 and 21, respectively), and thus are expected to harbo r little bias (< 20%) with respect 
to eccentricity according to the analysis of Shen Sz Turner (j2008l ). For this reason and the small 
size of the sample (2 systems) we neglect the eccentricity bias for SBls. 

One double-lined spectroscopic binary (SB2) has a measured orbit: PPl 15 (e = 0.42 it 0.05 
Basri &: Martinlll999l ). It was discovered as an SB2 in three epochs of RV data. The identification of 
SB2s in spectroscopic data is subject to different eccentricity selection effects than for SBls because 
changes in the velocity do not need to be measured to detect the two sets of absorption lines due 
to the binary. To quantify the eccentricity bias in detecting SB2s, we performed a Monte Carlo 
simulation similar to the one described above for visual binaries but with the detection criterion that 
the velocity shift between the binary components must be large enough to detect in high-resolution 
spectra at one of three random epochs. We fixed the period and semimajor axis to that of PPl 15 
(5.8 days, 0.03 AU)@ Figure [2] shows the resulting SB2 eccentricity bias predicted b y our simulation 
for a detection threshold of ARV = 11 km s~^ (i-e., 3x the spectral resolution of iBasri fc Martin 



19991 ). We find only a slight («3%) preference for the discovery of highly eccentric (e > 0.8) SB2s 
in such RV data. For comparison, we also show the results if the detection threshold were higher 
or lower by a factor of two, as would be the case for binaries with RV semiamplitudes different by 
a factor of two. Higher semiamplitude binaries are less biased, while for lower semiamplitudes the 
eccentricity bias becomes a larger effect (though still < 15%). Thus, it seems reasonable to neglect 
eccentricity bias for PPl 15, the one SB2 in our orbit sample. 

LSR J1610— 0040 (e = 0.444 it 0.017) is the o nly very low-mas s astrometric binary known to 
date, discovered by the USNO parallax progr am (IDahn et al.ll2008l ). It was later found to be an 
SBl in the multi-epoch RV data obtained by iBlake et al.l (120101 ). In principle, astrometric time 
series data is very similar to RV time series data and so would be subject to the same biases as 
for SBls. Thus, we assume that the same selection effects apply, namely that orbits with e > 0.85 
have an ~80% detection efficiency relative to less eccentric orbits, which are all equally accessible 
by astrometric monitoring. 

Finally, just one very low-mass eclips ing binary is known, 2MASS J0535— 0546 (e = 0.323 ± 
0.006), discovered bv IStassun et al.l (120061 ) in a photometric variability survey. In a study of ex- 
oplanet transit surveys. iBurkd (120081 ) showed that photometric transit/eclipse detections are only 



■^Note that our simulation on the detection of SB2s is nominally independent of an assumed mass ratio, as ARV 
does not depend on this quantity (ARV oc a sin iP^^Vl — ). Of course, very low mass ratio systems of all 
eccentricities are equally selected against as the secondary's absorption lines are more difficult to detect in the 
integrated-light spectrum. 
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slightly biased (< 10%) in eccentricity, with the exact correction depending on the underlying e 
distribution. Thus, we neglect any eccentricity bias for 2MASS J0535— 0546. 

Table H] lists the orbital periods, eccentricities, and masses for these five unresolved binaries' 
orbits. Combined with the visual binary sample, our sample comprises a total of 16 very low-mass 
binary orbits spanning orbital periods from 5.8 days to 8400 days. 



4. Discussion 



The eccentricity distribution and orbital period-eccentricity diagram for our sample is shown in 
Figure [3l Our sample shows a strong preference for modest eccentricities, with a median eccentricity 
of 0.34 (mean also 0.34) for the full sample and a median of 0.33 (mean of 0.33) for the visual 
binaries alone. The visual binary orbits show a broad range of eccentricities, from near-circular 
up to e = 0.830. At the shortest periods (P < 1000 days), we find an absence of eccentric orbits 
(e > 0.5), although this could be due to small number statistics. Thus, we find no strong evidence 
that the eccentricity distribution of very low-mass binaries changes over three orders of magnitude 
in orbital period (6 < P < 8 x 10'^ days). 

We have quantified the lack of a correlation between period and eccentricity by computing 
Spearman's rank correlation coefficient (rg) for our sample. We use a Monte Carlo approach, 
accounting for the period and eccentricity errors by drawing the data in each trial from the observed 
distributions (our MCMC chains when possible, Gaussians otherwise). The resulting mean and rms 
of the correlation coefficient was rs = —0.02 it 0.07, which is well below the critical value of 0.412 
needed to show correlation for our sample even at 90% confidence. The subset of visual binaries 
also shows no significant P-e correlation [rs = 0. 02 =fc 0.08, critical value 0.497 at 90% confidence) 
or a-e correlation [rs = 0.14 =b 0.09). In contrast. iKonopackv et al.l (j20ld ) reported a trend in e as 
a function of period for their sample of fifteen very low-mass visual binaries, although they did not 
quantify the significance of the correlation. We note that their observed trend is largely dependent 
on the two longest period orbits in their sample {P = 320 it 240 yr and 2000l^^ggQ yr) for which they 
found large eccentricities with large errors (e = 0.71 it 0.18 and 0.85^g 45! respectively). Both of 
these binaries were excluded from our sample by our criteria for minimizing eccentricity bias, and 
we also note that with only 2.2% and 0.4% observational phase coverage these orbit determinations 
an d their uncertainties are not likely to be robust (e.g., see the cautionary example of DF Tau AB 
inlSchaefer et al.lboodl. 



We have also investigated whether mass ratio is correlated with eccentricity for our sample. In 
Table m we list the best available mass ratio {q = M2/M1) estimates or measurement s for our sample 
binaries . For the yisual binaries, q is estimated from Lyon evolutionary models (IChabrier et al 



2000; iBaraffe et al.l 120031 ) given the measured total mass and individual luminosities using the 



method described in Section 5.4 of iDupuy et al.l (|2009bl ). This has been done in a consistent 
fashion for the sample presented here, so even if there are systematic errors in the evolutionary 
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models the relative ordering of the mass ratios should be relatively unaffected0 Fo r the eclipsing 



binar y 2MASS J0535— 0546, we used its directly measured mass ratio (0.64 it 0.04; IStassun et al 



20061 ). None of the three SBls in our sample has a well-determined mass ratio, so we excluded them 



from our mass ratio analysis. We also excluded one visual binary (LHS 2397aA B) because it has a 



very large uncertainty in its mass ratio estimate (0.71^q ]^^: lDupuv et al.ll2009d ) due to the fact that 
it is composed of a main-sequence star and a field brown dwarf. We computed Spearman's rank 
correlation coefficient in a Monte Carlo fashion as described above and found rs = —0.09 it 0.10, 
showing that there is no evidence for a correlation of eccentricity with mass ratio over the range of 
our sample (0.6 ^ ^ 1.0). In the absence of such a correlation, any mass ratio selection effects 
present in our sample should therefore not induce a systematic eccentricity bias. 



4.1. Comparison to Solar-type Binaries 



Duquennov &: Mavoij (|l99ll ) have examined in detail the eccentricity distribution of a well- 
defined sample of solar-type binaries. In the following analysis, we exclude from c onsideration three 
of the binaries in their sample that have an orbit quality grade of 4 or 5 assigned by lWorley &: Heintz 
(| 19831 ). which indicates an indeterminate orbit. Not surprisingly, t hese binaries (HP 16895AB , 
HD 78154AB, and HD 146361AB) are the longest pe riod orbits in thelDuauennov &: Mayor! (Il99ll ) 
sample, with orbital periods > 1000 yr. (We note that Raghavan et al. ( 20ld ) have recently updated 
this volume- limited, solar-type sample, but we cannot use this for comparison since they did not 
publish the eccentricity values of their sample binaries.) 



We find that the overall eccentricity distribution of the iDuquennoy &: Mayod (Il99ll ) solar- type 
binary sample is very similar to our orbit sample. To assess this quantitatively, we performed 
Kolmogorov-Smirnov (KS) tests in a Monte Carlo fashion. We accounted for eccentricity errors in 
our sample by drawing the data in each trial from the observed distributions (our MCMC chains 
when possible, Gaussians otherwise). We computed the KS statistic (D) and the corresponding 
probability for every trial. The resulting rms scatter in D was relatively small(< 0. 05)1 For 
the final probability that two samples were drawn from the same parent distribution, we used the 
mean of the probabilities from individual trials. Table [5] shows the res ulting probabilities compariii g 
subsamples of our orbits to subsamples of the solar- type binaries from lDuquennov fc: Mavorl (|l99ll ). 
In general, we find good agreement between the two samples' eccentricity distributions. However, 
a more detailed examination of these samples shows some significant differences. 



We note that a few visual binaries in our sample have mass ratios reported bv lKonopackv et al. I (|2O10l ) based on 
resolved radial velocities. However, they have unacceptably large uncertainties, with la ranges typically spanning a 
factor of > 5 in mass ratio except for Gl 569Bab {q = O.Tlly 131 i-e., 0.58-0.90 at la), which is still very large for our 
purposes. Thus, we have used our model-based q estimates for these binaries. 



^The definition of D is the maximum difference between the two cumulative distribution functions that are being 
tested. 
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Unlike the very low-mass binaries, solar- type binaries follow a strong correlation with orbital 
period such that short-period (< 10^ days) bina ries span a different range o f eccentricities than 
do long-period (> 10^ days) binaries (Figure H]). iDuquennov &: Mavorl (jl99ll ) suggested that this 
implies a different formation mechanism and/or dynamical history for these two period regimes 
of solar-t ype bina.ries. The same trend has also been reported for pre-main-sequence, solar-type 
binaries (jMathieulll994l . and references therein), implying that this period-eccentricity correlation 
is dete rmined early in a binary's fo rmation history. We confirm the significance of the correlation 
in the IDuquennov &: Mavoii (|l99ll ) sample by computing Spearman's rank correlation coefficient 
for solar-type binaries (excluding those that lie below the critical period for tidal circularization; 
-fcirc = 11-6 days). The resulting value of rs = 0.32 for this sample of 46 solar-type binaries 
gives a 98.6% confidence that P and e are correlated. As discussed above, our sample shows no 
such correlation, which is an indication that the solar and very low-mass populations are in fact 
dynamically distinct from one another. 



W e also note that among the large sample of solar-type binaries from iDuquennoy &: Mayor 
( 1991 ) there is not a single orbit with e < 0.1 beyond an orbital period of 22 days Ifl while there are at 
least 3 such orbits in our sample of 16 binaries (one more, 2MASS J1534— 2952AB, could also have 
e < 0.1, but the error in its eccentricity makes this uncert ain). This feat ure of solar- type systems 
has also been pointed out for pre-main-sequence binaries (jMathieulll994l . and references therein). 
We used a Fisher's exact te st to compute the pr o babili ty that this difference in e < 0.1 orbits 
between our sample and the IDuquennov &: Mavorl (jl99ll ) sample arises solely by chance (i.e., the 
P- value). In order to exclude any solar- type binaries that may be affected by tidal circularization 
but have not reached e = yet, we consider only binaries with P > 100 days (a less conservative, 
shorter period cutoff could be used but would not change the results). At these orbital periods, 
of 39 solar-type binaries have e < 0.1, while 3 of 14 very low-mass binaries have eccentricities 
confidently < 0.1. Comparing these two samples. Fisher's exact test finds that this difference is 
significant with a P- value of 0.0155 {2.2a). The P- value remains significant (< 0.05) for any choice 
of the eccentricity division point ranging from 0.10 to 0.15. This strong preference for near-circular 
orbits among very low-mass binaries provides additional evidence that their dynamical history may 
be different from solar-type binaries. 

We next consider the location in the P-e diagram of very low-mass binaries that are members 
of hierar chical triple systems a s com pare d to their solar-type c ounterparts. The interpretation 
given by IDuquennov Sz Mavoii (jl99ll ) and iRaghavan et al.l (|20ld ) is that solar-type binaries that 
are members of hierarchical triple or quadruple systems tend to lie near the upper envelope of 



®Solar-type binaries with P > 2 2 days and e < 0.1 ar e indee d known to exist, but none are inclu ded in the well- 
defin ed, volume-limited sample of iDuquennov fc Mavon (119911 ). For example, in the SB9 catalog (jPourbaix et al 



2004 ) 20% of orbits with P > 100 days and quality grades > 3 have e < 0.1. However, such spectroscopic binary 



catalogs represent a compilation of heterogeneous data not intended to be used for statistical studi es, as they likely 
harbo r significant selection biases (e.g., as discussed in Section [3.2p . In addition, we note that the iRaehavan et al 
( 2O10l ) solar-type sample contains 3 such binaries according to their Figure 14. 
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eccentricity at a given orbital p e r iod . However , exami nation of the P-e diagram of solar-type 
binaries presented in Figure 14 of lRaghavan et al.l ([20101) shows that such members of hierarchical 
systems actually occupy a wide range of eccentricities at periods > 100 days. In fact, members of 
hierarchical triple and quadruple systems are the only near-circular (e < 0.1) solar-type binaries 
at these periods, suggesting that their higher order nature may be the cause of their unusually low 
eccentricities. This is consistent wit h the fact that t he only such long-period, near-circular binary 
in the pre -main-sequence sam ple of iMathieul ()1994l ) is now known to have a tertiary companion 
(GW Ori; iBerger et al.ll201lh . In contrast to solar- type binaries, the higher order multiples in 
our sample do not populate the upper eccentricity envelope at any period. In fact, the highest 
eccentricity orbit in our sample is a true binary, and the lowest eccentricity orbit belongs to a 
member of a triple system. This latter datum may be partially consistent with the fact that near- 
circular (e < 0.1) solar-type binaries tend to be members of higher order multiples. However, there 
remains the mystery as to why such near-circular orbits are much more common among the true 
binaries in our sample than in the solar-type binary sample. 

Finally, we note that the shortest period binaries in our sample are not circularized, despite 
having shorter period s than the empirical li mit from field solar-type binaries (Pcirc = 11-6 days). 
This was discussed bv lBasri &: Martini (119991). who u sed brown dwarf structural models along with 
the circularization theory of IZahn &: Bouchetl (|l989l ) to show that PPl 15 {P = 5.8 days) is not 
expected to ever tidally circularize as its critical period is ~4.5 days. By extension, the orbit of 
the eclipsing binary 2MASS J0535— 0546 {P = 9.8 days), which has even lower masses and thus a 
shorter critical period, is also pred icted to never circul arize, so its nonzero eccentricity is not merely 
a symptom of its youth (~1 Myr; IStassun et al.ll2006l ). 



4.2. Comparison to Theory 



AmbartsumianI (|l937l ) showed that if binary companions are distributed in phase space solely 
as a function of energy, whether according to the Boltzmann distribution or any other arbitrary 
function, the cumulative distribution of eccentricities must be proportional to and the mean value 
of e should be equal to |. He argued that the observational data for wide stellar binaries at the time 
was in strong disagreement with this prediction, implying that binari es neither originate from no r 
relax via d ynamical interactions into a state of statistical equilibrium. iDuquennov &: Mavoii (|l99ll ) 
tested the lAmbartsumianI (|l937l ) distribu tion with modern data f or sol ar-type binaries, applying 
a correction due to eccentricity bias from [Harrington &: MiranianI (jl977l ). and claimed that the 
distribution described the wider binaries {P > 10^ days) reasonably well. However, their data 
still possess a severe lack of high-eccentricity binaries, especially after excluding the three binaries 
with >1000-yr orbits, as these happen to include two of the most eccentric orbits in their sample 
(e > 0.8). 

Our sample seems to also contravene the distribution, as high-e binaries are notably scarce. 
To quantify this, we performed KS tests as described in Section 14.11 Following the discussion in 
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Section 13. H we corrected for visual binary orbit bias by adding simulated data points uniformly 
distributed in each of the last two bins (e = 0.6-0.8 and 0.8-1.0). These bins contain and 1 
visual binary, respectively, and after accounting for the bias the number in each bin was and 2.4 
binaries, averaged over 10^ trials. The resulting probability that the null hypothesis is correct (i.e., 
the observed sample is drawn from the predicted distribution) was 9.1 x 10^^ for the entire sample 
and 2.2 x 10"'^ when considering only the visual binaries in our sample. Because of the potential 
impact of our assumed bias correction on the results, we also performed KS tests considering only 
eccentricities below 0.6, for which there is expected to be negligible observational bias. These tests 
gave probabilities of 5.3 x 10~^ and 7.1 x 10~^ for the full sample and visual binari es, respectively. 
Thus, we conclude that our sample of very low-mass binaries does not follow an lAmbartsumian 
(|l937l ) distribution over orbital periods of 5.8 to 8 x 10'^ days. 



We have also compared our sample to the p redict i ons of numerical simulations that strive to 
describe the formation of such low-mass binaries. iBatd (l2009al ) simulated the collapse of a 500 Mq 
molecular cloud, producing sufficiently la rge number s of stars and brown dwarfs (> 10'^) to do 
statistical comparisons with observations. iBatd (j2010l ) reported the parameters of sixteen binaries 
with total masses of < 0.2 Mq that were generated in this simulation, and we plot their orbital 
periods and ecce ntricit i es in F igure Our observed eccentricity distribution agrees very well with 
the simulation of lBatd (|2009al ). with a 83% probability of both samples being drawn from the same 
parent distribution according to a KS test (92% for visual binaries alone), and this result is not 
changed if we consider only those sim ulated binari es within our sample's semimajor axis upper 
limit of 3.7 AU. However, we note that iBatd (l2009al ) does not produce any high-e (> 0.6) binaries 
over the same orbital period range as we observe, as all of his high-e binaries have very long orbital 
periods (> 10^ days). This could be due to the small number statistics of the simulation results, but 
if not then it is suggestive of a disagreement with the observations that is only made stronger by the 
fact that the orbit bias correction implies a total of ~3 binaries with e > 0.6 for our sample. This 
discrepancy could be due to the fact that at the end of the simulation the iBatd (|2009al ) binaries are 
still evolving toward shorter orbital periods by dynamical interactions with other protostars and 
their own accretion disks. Thus, the high-e systems could ultimately reach the orbital periods of 
our sample. However, their eccentricities would also likely damp down in the process, and thus the 
discrepancy with observations would remain . Another sig nificant concern is that the very low-mass 
stars and brown dwarfs in the simulation of iBatd (l2009al ) are for i ned lar gely by being ejected from 
their initial gas reservoirs through dynamical interactions. iBatd (j2009bl ) showed that by including 
a more realistic equation of state and radiative feedback, such ejection ev ents become much less 
common, which resolves the problem that previous simulations (including iBatd l2009al ) produced 
far too many brown dwarfs. It remains to be seen how this improvement in the input physics will 
impact the orbital elements of binaries produced in these simulations. 



^Note that these are results from the "rerun" simulation, whi ch had a smal ler sink radius of 0.5-AU, no gravitational 
softening, and a shorter run time than the main calculation in iBatd (|2009al V 
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In recent simulations bv lStamatellos &: WhitworthI ((20091), very low-mass stars, brown dwarfs, 
and planetary-mass objects are formed via gravitational instability in the outer parts (> 100 AU) of 
massive circumstellar disks (Mdisk/-^* = 1-0). The orbital periods and eccentricities of the twelve 
binaries with Mtot < 0.2 that were generated in these simulations are shown in Figure HI These 
binaries are a mix of three field binaries (those that are ejected from the disk), and nine binaries 
that remain bound to the star in a hierarchical triple. The eccentricities of these two subsets are 
not statistically different from each other as shown by a KS test {D = 0.5, 38% chance of being 
drawn from the same parent distribution), and so we simply compare the full set of simulated 
binaries to our sample. The simulated binaries have significantly higher eccentricities than we 
observe, with our KS test giving a probability of only 9.9 x 10~^ that the two samples were drawn 
from the same parent distribution. This result is not changed significantly if we only consider 
simulated binaries within our s ample's semimajor axis upper limit of 3.7 AU. As discussed by 



Stamatellos &: WhitworthI (j2009l ). the binary components formed in these simulations have their 



own disks, and after the initial 20 kyr their simulations ignore the gas in order to follow the 
long term evolution of the system. The tidal interactions associated with these disks that are not 
account ed for in the s imulations at late times could damp the eccentricities of the simulated binaries. 
Indeed, Bate (j2009al ) found that including such dissipative interactions was critical in reproducing 
the stellar eccentricity distribution. Thus, including such additional in put physics may solve the 
disagr eement between our observed eccentricities and the simulations of IStamatellos Whitworth 
([20091). On the other hand, such large eccentricities may simply be a hallmark of forming in the 
frenetic dynamical environment of a fragmenting circumstellar disk. 



Conclusions 



We have considered for the first time the eccentricity distribution of a large sample of orbits 
for very low-mass star and brown dwarf binaries. We have examined potential observational biases 
in this sample, and for visual binaries we show through Monte Carlo simulations that by choosing 
appropriate selection criteria all eccentricities are equally represented (< 5% difference between 
input and output eccentricity distributions). We find that our sample of very low-mass binaries 
populate a broad range in eccentricity from nearly circular to highly eccentric (0.03 < e < 0.83), 
showing a strong preference for modest eccentricities (the mean and median of the full sample is 
e = 0.34). The binaries among our sample that are members of hierarchical triple systems show no 
special preference for their eccentricity as compared to the other binaries, consistent with a common 
formation mechanism for both po pulations. The overall eccent ricity distribution of our sample is 
similar to solar-type binaries from iDuquennoy &: Mayod (|l99ll ). However, very low-mass binaries 
lack any correlation between eccentricity and period or semimajor axis, which is one of the key 
features of the solar-type population. In addition, near-circular orbits appear to be quite common 
among very low-mass binaries, as at least 3 of our 16 sam ple binaries have e < 0.1. H owever, such 
orbits are apparently rare among solar-type binaries (in the lDuquennoy &: Mayoiill99ll sample, none 
of the 42 binaries beyond an orbital period of 22 days have e < 0.1). This difference between the 
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two samples is statistically significant according to Fisher's exact test (P- value of 0.0155). Thus, 
the orbital data currently available are suggestive of a different formation mechanism or dynamical 
history for very low-mass binaries as compared to solar- type stars. 

We have also tested formation models by comparing our sample of very l o w-mas s binaries to 
those produced in numerical simulations. The cluster formation model of iBatd (|2009al ) predicts an 
overall eccentricity distribution that is well-matched to the data. However, there is a potential dis- 
crepancy with the observations due to the lack of any eccentric binaries (e > 0.6) at orbital periods 
corresponding to our sample. This may be caused by the same physical mechanism that is needed 
to bring binaries formed in the simulation to these orbital periods, namel y tidal damping by the ac- 
compa nying gas disks. In contrast, the gravitational instability model of lStamatellos &: Whitworth 
([20091) produces too many high-e binaries, with an absence of any modest eccentricity (e < 0.3) 
orbits. This may be due to th e fact that tidal dissipation by gas disks is not included in the 
Stamatellos &: WhitworthI (|2009l ) simulation, causing the resulting eccentricity distribution to be 
inconsistent with our sample at high significance. A third possibility is that the observed discrep- 
ancies are simply due to multiple formation processes being responsible for producing the field 
population of very low-mass binaries, so no one model matches all the data. A further complica- 
tion is that binaries formed in a cluster environment are likely to undergo additional dynamical 
processing after their initial eccentricities are set, and thus A^-body simulations are needed to assess 
the impact of such effects on the final eccentricity distribution. 

To improve on our observational tests, future efforts should aim to cover a wider range of orbital 
parameter space, particularly expanding the sample at short orbital periods {P < 10^ days). The 
discovery of such tight binaries is only currently feasible from RV monitoring or the "spectral 
binary" technique (i.e., identifying the signature o f a fainter companion o f disparate spectral type 
in the integrated light spectrum). For example, iBurgasser et al.l (|2008l ) used this t echnique to 



identify 2MASS J0320— 0446 as a binary independent of its identification as an SBl (jBlake et al 



20081 ). However, even for such spectral binaries the determination of orbits requires either precise 
RV measurements or resolved astrometry. This will require new capabilities (precision near-infrared 
radial velocities for L and T dwarfs) as well as new facilities (e.g., TMT/GMT/E-ELT and JWST to 
resolve the much smaller semimajor axes). Fortunately, future binaries that will fill in the bottom 
two orders of magnitude in orbital period will not require the same duration of patient orbital 
monitoring that has been necessary over the last decade to develop the current sample of ~3-30-yr 
period visual orbits. 

We have taken the first look at what clues the ensemble of orbital eccentricities of very low- 
mass stars and brown dwarfs reveal about their formation. This work follows on previous statistical 
studies of the visual binary population that found distributions of mass ratio and semimajor axis 



that were significantly different than for solar-type binaries (e.g.. lClose et al.ll2003l : IBurgasser et al 



20071 ) . The theoretical implication that disks may be critical in shaping the eccentricity distribution 
is in accord with the fact that eccentricity is closely associated with a binary's angular momentum. 
This potential relationship is intriguing, as it suggests the possibility for disk-related observational 
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tests. For example, if high-eccentricity binaries are simply those that lost their disks early in their 
formation, then perhaps the binary components rotate more rapidly than lower eccentricity binary 
components due to a lack of disk-braking. Observational evidence has also shown that the forma- 
tion environment may play an important role in determining the orbits of binaries, with dynamical 
interactions in the dense Orion Nebula Cluster likely responsible for sculpting the semimajor axis 
distribution of stellar binaries (jReipurth et alJ 120071 ) while binaries in the sparse r Taurus region 
seem to be able to survive at much wider separations (iKraus &: Hillenbrand! l2009l ) . Further theo- 
retical modeling of the binary formation process taking such effects into account will enable us to 
better understand the physical processes influencing observed orbital eccentricities. 



Our research has employed NASA's Astrophysical Data System; and the SIMBAD database 
operated at CDS, Strasbourg, France. We are grateful to Matthew Bate for providing us with the 
orbital elements of the lowest mass binaries in his simulation. We have benefitted from discussions 
with Adam Kraus about binary formation, John Johnson about RV selection effects, and Tsevi 
Mazeh about orbital dynamics. TJD and MCL acknowledge support for this work from NSF 
grants AST-0507833 and AST-0909222. TJD acknowledges support from Hubble Fellowship grant 
HST-HF-51271.01-A awarded by the Space Telescope Science Institute, which is operated by AURA 
for NASA, under contract NAS 5-26555. Finally, the authors wish to recognize and acknowledge 
the very significant cultural role and reverence that the summit of Mauna Kea has always had 
within the indigenous Hawaiian community. We are most fortunate to have the opportunity to 
conduct observations from this mountain. 

Facilities: Keck H Telescope (LGS AO, NGS AO, NIRC2) 



A. Converting Projected Separation to Semimajor Axis 



A general problem in the study of visual binary stars is that orbit determinations are often 
infeasible, so the true semimajor axis (a) of the binary must be estimated from the projected 
separation on the sky (p). In the usual approach, p is multiplied by a correction factor that 
accounts for the fact that the observed separation will typically be smaller than the true semimajor 
axis. This is always true for circular orbit s, but eccent r ic orb its can be observed at much wider 
separations than the actual semimajor axis. lAbt &: Levvl (1197611 used a correction fact or of a/p = ^, 
computed analytically by assuming edge-on circular orbits. iFischer Marcyl (jl992l ) performed a 
Monte Carlo simulation that assumed random inclinations and reported 1.26 as their time-averaged 
correction factor. Unfortunately, they did not describe the input eccentr icity distribution used nor 
did they report confidence limits for the correction factor. I Torres! (jl999l ) performed several Monte 
Carlo simulations assuming random viewing angles and a variety input eccentricity distributions. 
He showed that the conversion factor not only varies greatly for different eccentricities, but that 
the probability distributions for these values are quite broad and can display one or more peaks. 
Thus, it is important to account for these large uncertainties when converting p to a. 
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With empirical eccentricity distributions in hand for solar-type and very low-mass binaries, we 
have undertaken Monte Carlo simulations to compute appropriate conversion factors for binaries 
drawn from these populations. We assumed uniformly distributed observation times {0 < t — to < 
P), uniform arguments of periastron (0° < ui < 360°), and randomly oriented orbits {p{i) oc sini). 
From 10^ Monte Carlo trials we computed observed projected separations (p) and thus arrived at 
a distribution of conversion factors (a/p)- The results for various assumptions about the input 
eccentricity distribution are given in Table [6l In Figure [5] we show the probability distributions of 
these conversion factors. For comparison, we also show the resulting distribution if e is held constant 
in the simulations. In this case, the distribution has two peaks: one at 1/(1 -|- e) corresponding to 
its apastron passage and one at 1/(1 — e) corresponding to periastron. The empirical distributions 
of a/p have many small localized peaks due to the fact that they are composed of many individual 
eccentricity measurements that each contribute two peaks. The effect is more pronounced for the 
smaller sample of the thirteen very low-mass binaries as compared to the 46 solar- type binaries. 

We have further considered the impact of discovery bias (see Section [3T]l on these distributions. 
The resolution limit of an imaging survey will preclude the discovery of binaries when they are at 
very small separations, which results in truncating the tail of conversion factors at large a/p. In 
Table [6] we give the results when considering two cases of discovery bias. In the moderate case, 
the inner working angle (IWA) is a factor of two smaller than the true semimajor axis (IWA 
= a/2), which truncates the conversion factor distribution at a/p = 2. In the severe case the 
resolution limit is equal to the semimajor axis, so only eccentric binaries can be discovered, and the 
conversion factors are never greater than a/p = 1. These truncation values are shown as vertical 
lines in Figure [H 
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Fig. 1. — Top: Fraction of visual binaries detected in Monte Carlo simulations of a single-epoch 
imaging survey, where the input eccentricity distribution is uniform. Different colored lines show 
simulations in which the semimajor axis (a) was fixed at different multiples of the resolution limit, 
or inner working angle (IWA). When the semimajor axis lies inside the IWA highly eccentric orbits 
are strongly preferred, and as the IWA goes to zero the eccentricity bias disappears. Bottom: 
Eccentricity bias for binaries with late-M (solid) and L dwarf (dashed) primaries as determined 
from Monte Carlo s imulations a ssuming a magnitude- limited sample with semimajor axis drawn 
from the log- normal I AllenI (|2007l ) distribution. The resolution limit for both simulations was chosen 
to be the worst-case from among all surveys that discovered the binaries in our sample (0'.'12 for 
late-M dwarfs and O'.'OT for L dwarfs). The resulting worst-case eccentricity bias is modest, only 
slightly (~5%) favoring the discovery of highly eccentric (e > 0.8) and nearly circular (e < 0.2) 
binaries. 
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Fig. 2. — To'p: Fraction of SB2s detected in Monte Carlo simulations of a three-epoch RV survey, 
where the input eccentricity distribution is uniform. The solid black line is a simulation designed to 
imitate the case of PPl 15, the only SB2 in our orbit sample. Dotted colored lines show simulations 
in which the detection threshold (ARV = 11 km s~^) was changed by a factor of two. The RV 
observations that discovered PPl 15 only slightly (~3%) favored the detection of highly eccentric 
(e > 0.8) binaries. This effect would be larger (~15%) for binaries with a factor of two lower RV 
semiamplitudes. Bottom: Same as top panel, except that the dotted colored lines show variations 
in the eccentricity bias due to having fewer than the three epochs of RV data that were available 
for the detection of PPl 15. Surveys with only one or two RV epochs are less sensitive to SB2s at 
all eccentricities, and they are more biased toward the discovery of highly eccentric binaries. 
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Fig. 3. — Top: Eccentricity distributions for our sample of eleven very low-mass visual binaries 
(left) and the full sample including spectroscopic binaries (right). The dotted histograms show the 
orbit-fitting bias correction we applied to visual binaries in the last two bins (e > 0.6). Histograms 
were generated from 10'* randomly drawn eccentricity values for each binary corresponding to the 
measurement error, so some bins may have a non-integer value. Bottom: Eccentricity plotted as a 
function of semimajor axis for the visual binary sample (left) and as a function of orbital period for 
the full sample including spectroscopic binaries (right). Triangles denote members of hierarchical 
triple systems, and circles denote true binaries. 
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Fig. 4. — Eccentricity distributions (left) and orbital period -eccentricity diagrams (right , open 
symbols) for various samples. Top: Solar- type binaries from iDuquennov &: Mavoii (|l99ll ). with 
three low-quality, >1000-yr orbits excluded. Tidally circularized orbits (P < Pcirc) are plotted 
separately in the histogram (hatched region). Triangles and squares denote members of hierar- 
chical triple and quadrup le systems, respectively. Middle: Cluster formation model simulated 
binaries fro mlBa.tel (l2009al ): eccentricities for this sample of very low- mass simulated binaries are 
reported by Batel ( 201C ). B ottom : Circumstellar disk fragmentation model simulated binaries from 



Stamatellos &: Whit wort hi (120091 ). with triangles denoting binaries that remain bound to the sim- 



ulated host star. The very low-mass binaries in our sample are shown as filled red symbols on all 
P-e diagrams for comparison. 
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Table 1. New Keck AO Astrometry 



Epoch (UT) p (mas) PA (°) Filter Am (mag) 



LP 415-20AB (LGS) 



2008 Jan 15 89.3 ± 0.4 56.0 ± 0.5 Kg 0.60 ± 0.07 

2008 Sep 8 107.5 ± 0.8 61.0 ± 0.3 Ks 0.589 ± 0.013 

2009 Sep 28 129.81 ±0.15 66.70 ± 0.06 Ks 0.524 ± 0.022 

2010 Jan 9 134.96 ±0.19 68.05 ± 0.05 Kg 0.548 ± 0.016 



HD 130948BC (NGS) 



2010 Jan 9 89.0 ± 0.4 318.21 ± 0.25 Kg 0.212 ± 0.023 
2010 Mar 22 100.15 ± 0.15 316.95 ± 0.07 K 0.206 ± 0.007 



-28- 



Table 2. Updated Orbit Determinations 



Parameter LP 415-20AB HD 130948BC 



P (yr) 


14.4 ±0.4 


10.28 ±0.17 


a (mas) 


93.5+1^ 


124.5 ± 1.4 


e 


0.7081°:°}^ 


0.176 ± 0.006 


in 


59.9 ±2.2 


95.7 ±0.2 


nn 


261.511^ 


313.39 ±0.12 


un 


351 ±7 


244 ±3 


To (MJD) 


53505 ± 80 


54613 ± 22 


Mtot" (Mo) 




0.1095 ±0.0022 


(DOF) 


6.83 (7) 


9.38 (13) 



^The total mass accounts for the uncertainty in 
the parallax. LP 415-20AB has no measured par- 
allax and thus no direct mass measurement. 



Table 3. Visual Binaries with Published Orbits 



Name 




Sp. Type 


Distance 






Orbit 




Orbit Quality 


IWA 


IWA/a 






Primary 


d (pc) 


Ref. 


a (AU) 


e 


Rcf. 


Atobs/P 


cjm/M 


(") 




Gl 569Bab 




M8.5 


9.65 ±0.16 


16 


0.923 ±0.018 


0.316 ±0.005 


8 


334% 


— 3% -1-4% 


0'.'05 


0.52 


2MASS J09204 


-3517AB 


TO 


25.41°;^ 


9 


^^■'"-0.05 


0.217 ±0.019 


9 


144% 


3% 


0'.'06 


0.91 


LP 349-25AB 




M8 


13.2 ±0.3 


10 


1.94 ±0.04 


0.051 ±0.003 


8 


76% 


0.6% 


O^'IO 


0.68 


DENIS J2252- 


1730AB 


L6 


15.0 ±0.3 


9 


rj + O. 5 

^•^-0.2 


„ LA 1 1 

'-'•^^-0.08 


9 


34%-56% 


5% 


0'.'07 


0.48 


HD 130948BC 




L4 


18.17 ±0.11 


16 


2.26 ±0.03 


0.176 lb 0.006 


1 


73% 


0.8% 


0'.'06 


0.48 


2MASS J2132^ 


-1341AB 


L5 


27.5 ± 0.5 


9 


2.24 ± 0.10 


'-'•'-'^-0.02 


9 


35% 


5% 


0'.'05 


0.61 


e Indi Bab 




Tl 


3.624 ± 0.004 


16 


2.43 ± 0.11 


0.53 lb 0.03 


2 


38% 


0.9% 


0'.'08 


0.12 


LP 415-20AB 




M7.5 


[27 ± 4]'^ 


1 


[2.5 ± 0.4] 


„ _.„,-,+0 016 


1 


55% 


-7%, +9% 


0'.'12 


1.28 


2MASS J07464 


-2000AB 


LO 


12.21 ± 0.04 


5 


2.899 ± 0.013 


0.487 ± 0.004 


1,11 


68% 


2% 


0'.'06 


0.25 


LHS 2397aAB 




M8 


14.3 ± 0.4 


13 


3.09 ± 0.10 


0.350 ±0.005 


7 


83% 


3% 


0'.'12 


0.56 


2MASS J1534- 


-2952AB 


T5 


13.59 ±0.22 


15 


Q 9+0.4 
■""•^-O.S 


in+0.09 

'J- '-'-0.06 


1,11 


32%-44% 


2% 


0'.'06 


0.26 


LHS 1070BC 




M8.5 


7.72 ±0.15 


3 


3.56 ±0.07 


„ „oo+0.003 

U.UOO_o QQ2 


1,14 


82% 


0.9% 


O'.'IO 


0.22 


LHS 1901AB 




M6.5 


12.9 ±0.5 


12 


3.70 ±0.16 


n 104+0.025 
u.iy4_o 021 


8 


39% 


1.4% 


O'.'IO 


0.35 


2MASS J21404 


-1625AB 


M8.5 


[26±5]=' 


1 


[3.7 ± 0.7] 


36+° " 


1,11 


38%-57% 


-27%, +53% 


0'.'12 


0.83 


2MASS J08504 


-1057AB 


L7 


38+^ 


17 




0.64 ±0.26 


11 


10%-49% 


100% 


0'.'06 


0.48 


2MASS J1728^ 


-3948AB 


L5 


24 1+2'^ 


17 


5.3 ±0.8 


o.28l°:i 


11 


20%-47% 


-27%, +167% 


0'.'06 


0.27 


2MASS J2206- 


-2047AB 


M8 


26.7t2;i 


4 


c O+0.8 
0-°-0.7 


0.25 ±0.08 


14 


20%-28% 


2.0% 


0'.'12 


0.56 


2MASS J18474 


-5522AB 


M6.5 


[25±5]=' 


1 


[5.8 ± 1.5] 




11 


9%-23% 


-72%, +194% 


0'.'08 


0.34 


2MASS J17504 


-4424AB 


M7.5 


[29±8]=' 


1 


[21 ± 13] 


0.71 ±0.18 


11 


1.3%-9% 


60% 


0'.'12 


0.16 


2MASS J1426+1557AB 


M8.5 


[25±5]=' 


1 


[60 ± 40] 


0.85«i; 


11 


0.2%-20% 


-100%, +73% 


0'.'12 


0.05 



Note. — This table includes all resolved binaries with a total mass < 0.2 Mq for which an orbit determination has been published. Orbit quality metrics are given for 
each binary: (1) fraction of orbital phase observed {/^t^hsl P)^ if the period uncertainty is > 10% the Icr range is given; and (2) fractional Itr error in the total mass from 
the orbit fit alone {um /M), i.e., independent of the uncertainty in the distance. High quality orbits have a large fraction of the orbit covered and/or a low mass error. 

^Spectrophotometric distance estimate derived following the method used in the Appendix of lLiu et al. I ll20ld'l . 



Refe rence s. — (l)This work; (2) ICardoso et~al] 1 20091 ') ; 
ll2009d'): (81|Duduv et aD ||201C); O*) IDuduv " | 201C ): (11) 
ri4) ISeifahrt et al.1 l|20o1 ): a5l iTinnev et al.1 \200± . (16) 



(3'llCosta et al.l ||2005|'): ( 4) ICosta et al.l ll2006f): (5')lDahn et al.l 1I2OO2I'): (6)IDupuv et al.l l l2009al'l; (7)IDupuy et al.l 
Konopackv et al.l jpOld ) ; (W) ICatewood fc CobanI l|2009l l: (12) iLepine et al.l 1120091 ): C13') iMonet et al.l |l993 h 
van LeeuwenI l|2007h : fl7l IVrba et al.l l|2004l l. 



Table 4. Very Low-Mass Binary Orbit Sample 



Name 


Int .-Light 


Binary 


p 




Mtot 


Mass ratio 


Orbit 




Sp. Type 


Type 


(days) 




(Mq) 


(o = M2/M1) 


Ref. 


PPl 15 


M6.5 


S2 


5.8 ±0.3 


0.42 ±0.05 


[0.124 ±0.004]'^ 


[0.85 ± 0.05] 


2 


2MASS J0535-0546 


M6.5 


E 


9.77962 ± 0.00004 


0.323 ± 0.006 


0.088 ± 0.005 


0.64 ± 0.04 


12 


2MASS J0320-0446 


MS 


SI 


246.9 ±0.5 


0.067 ±0.015 


[0.15 ± 0.03]=' 


[0.60-0.87] 


3 


LSR J1610-0040 


d/sdM7 


A, SI 


607 ±4 


0.444 ±0.017 


[0.162 ± O.OIS]'' 


[0.61-0.86] 


5 


Gl 569Bab 


M8.5 




864.5 ± 1.1 


0.316 ±0.005 


n 140+0-009 


[0.866ini!] 


7 


/-IT, „ TT „ O 

Una tia o 


M5.75 


SI 


1 r\r\r\ i i on 

1900 ± 130 


0.59 ± 0.22 


To 1 cr 1 ri ool a 

[0.15 ± 0.02]" 


[0.3-0.7] 


9 


LP 349-25AB 


M8 


V 


2834 ± 15 


0.051 ±0.003 




[0-872ini|] 


7 


HD 130948BC 


L4 




3760 ± 60 


0.176 ±0.006 


0.1095 ± 0.0022 


[0.948 ± 0.005] 


1 


€ Indi Bab 


T2.5 




4090 ± 180 


0.53 ±0.03 


0.116 ± 0.001 


[0.63 ± 0.02] 


4 


2MASS J2132+1341AB 


L6 


V 


40001^00 


34+0-03 


[0.13 ± 0.02]=^ 


[0.76 ± 0.03] 


8 


DENIS J2252-1730AB 


L7.5 


V 


40001^500 


Q -2+0-11 
'^■^^-0-08 


[0.12 ± 0.02]'' 


[0.70 ± 0.03] 


8 


2MASS J0746+2000AB 


L0.5 


V 


4640 ± 40 


0.487 ±0.004 


0.151 ± 0.003 


[0.94 ± 0.02] 


1,10 


LHS 2397aAB 


MS 


V 


5190 ±40 


0.350 ± 0.005 


0.146 ± 0.014 


[O.TllH!] 


6 


LHS 1901AB 


M7 


V 


5880 ± 180 


0.830 ± 0.005 


194+0-025 


[0.958ini!l 


7 


LHS 1070BC 


M8.5 




6220 ± 16 


033+0-003 

U.UOO_Q QQ2 


0.156 ± 0.009 


[0.84 ± 0.05] 


1,11 


2MASS J1534-2952AB 


T5 


V 


8400115°° 


10+00^ 
'^■^'^-0.06 


0.061 ± 0.003 


[0.936 ± 0.012] 


1,10 



CO 

o 



Note. — This table includes visual binaries from Table [3] that pass our eccentricity de-biasing criteria (0.9 AU < a < 3.7 AU and 
IWA/a < 0.75), as well as all unresolved binaries with a total mass < 0.2 Mq for which an orbit determination has been published. 
The type of each binary is given: astrometric (A), spectroscopic single-lined (SI) or double-lined (S2), eclipsing (E), or visual (V). 
All table values encl osed in brackets "[]" are estimates based on models. Mass ratio estimates for visual binaries are based on Lyon 
evolutionary models I Chabrier et al.lbood baraffe et al.ll2003l). while the estimates for the SB Is and SB2 are taken from the literature 
l lBasri fc Marti'nlll999l: IPahn et ahlbOOSl : iBurgasser fc Blakell2009l : Ijoergens et al.ll2010h . 

''No dir ect mass measurement is available, so es t imates from the literature based on evolutionary models and/or ancillary data 
are given iBasri fc Martini 119991 : ISiegler et al] |2003| : [Reid et al.l |2006| : ISiegler et al]|2007l: iDahn et al.ll2008l : iBurgasser fc Blakel ISooil : 
Ijoergens et al.ll2010l) . 

^Binary is the short-period component in a hierarchical triple system. 



Ref erenc es. — (1) This wo rk; (2) iBasri fc Martini lll999l'): (3) ISlake et al] ll2010h: (4) ICardoso et al.l ||2009|'): (51 IPahn et al.l 
| |200S|'I: (6) IDupuv et al] l l2009cll: (71 IDupuv et al.l toidi) : (8) lDupuvll2010l) : (9) Ijoergens et al.l boid) : (10) iKonopackv et al.l boT^ . 
(lll lSeifahrt et al.l 1 12008) : (12) Istassun et al.1 l l20od)" 
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Table 5. Results of KS Tests 



Very Low-Mass Binaries 

Comparison Sample All Visual only e < 0.6 





Prob. 


D 


Prob. 


D 


Prob. 


D 


Solar- type binaries (P > Pcirc) 


0.50 


0.23 ±0.04 


0.44 


0.28 ±0.04 


0.57 


0.23 ±0.04 


Solar-type binaries (Pcirc < P < 10"^ days) 


0.75 


0.22 ± 0.02 


0.82 


0.22 ±0.04 


0.74 


0.23 ±0.03 


Solar-type binaries (P > 10^ days) 


0.25 


0.31 ±0.05 


0.22 


0.36 ±0.04 


0.46 


0.28 ±0.03 


f Ambartsumian 1937) 


9.1 X lO"'^ 


0.55 ±0.04 


2.2 X 10^3 


0.52 ±0.04 


5.3 X lO"'' 


0.71 ±0.03 


Bate (2009a) 


0.83 


0.20 ± 0.04 


0.92 


0.19 ±0.03 


0.64 


0.26 ±0.04 


Bate (2009a. a < 3.7 AU) 


0.95 


0.22 ± 0.04 


0.97 


0.22 ±0.03 


0.91 


0.25 ± 0.05 


Stamatellos & Whitworth (2009) 


9.9 X 10-3 


0.59 ±0.05 


0.034 


0.56 ±0.05 


0.19 


0.52 ± 0.04 


Stamatellos & Whitworth (2009. a < 3.7 AU) 


0.028 


0.56 ±0.05 


0.064 


0.53 ±0.05 


0.43 


0.45 ± 0.05 



Note. — Each table entry shows the probability that our sample of very low-mass orbits was drawn from the same parent distribution 
as the comparison sample along with the D statistic and its rms scatter over 10* Monte Carlo trials. 
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Table 6. Conversion Factors from Projected Separation to Semimajor Axis 



Assumed e Distribution 


Correction factor 






Median 


68.3% c.l. 


95.4% c.l. 


No discovery bias 


Uniform (0 < e < 1) 


1.10 


0.75, 2.02 


0.57, 5.53 


Circular (e = 0) 


1.15 


1.01, 1.85 


1.00, 4.72 


e?" f Ambartsumian 1937) 


1.02 


0.67, 2.07 


0.55, 6.00 


Solar-type binaries (P > Pcirc) 


1.15 


0.81, 2.01 


0.64, 5.11 


Solar-type binaries (Pcirc < P < 10^ days) 


1.16 


0.83, 2.00 


0.66, 5.04 


Solar-type binaries (P > 10"^ days) 


1.14 


0.78, 2.02 


0.61, 5.29 


Very low-mass (all binaries) 


1.16 


0.84, 1.97 


0.66, 5.08 


Very low-mass (visual binaries) 


1.16 


0.85, 1.97 


0.66, 5.11 


Moderate discovery bias (IWA = 


a/2) 




Uniform (0 < e < 1) 


1.02 


0.72, 1.46 


0.57, 1.89 


Circular (e = 0) 


1.11 


1.01, 1.46 


1.00, 1.88 


f Ambartsumian 1937) 


0.91 


0.65, 1.42 


0.55, 1.88 


Solar-type binaries (P > Pcirc) 


1.06 


0.79, 1.51 


0.63, 1.90 


Solar-type binaries (Pcirc < P < 10^ days) 


1.07 


0.80, 1.51 


0.65, 1.90 


Solar-type binaries (P > 10^ days) 


1.04 


0.76, 1.50 


0.60, 1.90 


Very low-mass (all binaries) 


1.08 


0.82, 1.51 


0.66, 1.90 


Very low-mass (visual binaries) 


1.08 


0.82, 1.50 


0.64, 1.90 


Severe discovery bias (IWA = 


aY 




Uniform (0 < e < 1) 


0.79 


0.63, 0.94 


0.54, 0.99 


e?" (Ambartsumian 1937) 


0.74 


0.61, 0.91 


0.53, 0.99 


Solar-type binaries (P > Pcirc) 


0.83 


0.70, 0.95 


0.60, 0.99 


Solar-type binaries (Pcirc < P < 10^ days) 


0.84 


0.71, 0.95 


0.61, 0.99 


Solar-type binaries (P > 10"^ days) 


0.82 


0.68, 0.94 


0.58, 0.99 


Very low-mass (all binaries) 


0.84 


0.72, 0.96 


0.58, 0.99 


Very low-mass (visual binaries) 


0.85 


0.71, 0.96 


0.58, 0.99 



^When IWA > a, circular orbits are never detectable. 



